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Abstract. We study the problem of feedback control for skew-symmetric and skew- 
Hamiltonian transfer functions using skew-symmetric controllers. This extends work of 
Helmke, et al., who studied static symmetric feedback control of symmetric and Hamil- 
tonian linear systems. We identify spaces of linear systems with symmetry as natural 
subvarieties of the moduli space of rational curves in a Grassmannian, give necessary and 
sufficient conditions for pole placement by static skew-symmetric complex feedback, and 
use Schubert calculus for the orthogonal Grassmannian to count the number of feedback 
laws when there are finitely many. Finally, we also construct a skew-symmetric linear 
system and poles with only real feedback laws. 



1. Introduction 

Many fundamental questions about output feedback pole assignment for general linear 
systems have been answered by appealing to algebraic geometry, and more specifically to 
the geometry of Grassmann manifolds. This body of work has led to important contribu- 
tions in systems theory: Hermann and Martin gave necessary and sufficient conditions for 
complex static output feedback control [131 122l [21], Brockett and Byrnes used Schubert 
calculus to count the number of pole- assigning feedback laws [3], and then Rosenthal [2B] 
and Ravi, Rosenthal, and Wang [251 El] solved these problems for complex dynamic com- 
pensators using quantum Schubert calculus. For a description of the earlier literature, we 
recommend [2| • This line of work on complex feedback has led to a solution of the problem 
of pole-assignment in the real case, some of which is found in [71 [271 [32]. Likewise, it has 
influenced work in algebraic geometry [151 [291 [H] ) some of which is surveyed in [30] . 

The Lagrangian Grassmannian and orthogonal Grassmannian are subsets of the usual 
Grassmannian, and in principle they should also appear in systems theory. This was real- 
ized by Helmke, Rosenthal, and Wang [12], who studied the control of linear systems with 
symmetric and Hamiltonian state-space realizations by static symmetric output feedback. 
They gave necessary and sufficient conditions for pole placement by static symmetric 
complex feedback, linking this problem to the Schubert calculus on Lagrangian Grass- 
mannians, and then used this link to count the number of feedback laws when they are 
finite. 



1991 Mathematics Subject Classification. 93B55 (14M15 93B27). 

Key words and phrases, pole placement, feedback control, orthogonal Grassmannian, Lagrangian 
Grassmannian, skew-symmetric matrix. 

Hillar supported in part by an NSF Postdoctoral Fellowship and an NSA Young Investigator grant. 

Sottilc supported in part by the Institute for Mathematics and its Applications, Institut Mittag-Lcffler, 
and NSF grants DMS-0701050 and DMS-1001615. 

1 



2 



CHRISTOPHER J. HILLAR AND FRANK SOTTILE 



We continue this line of research. We first identify spaces of hnear systems with a natural 
symmetry as certain subvarieties of the space of rational curves in a Grassmannian. More 
specifically, we consider linear systems with McMillan degree n whose transfer function 
G{s) (a square matrix of rational functions) has one of the following four symmetries: 

(1) G{sY = G{s) symmetric, 

(2) G{s)^ = G{—s) Hamiltonian {n must be even), 

(3) G{s)^ = — G(— s) skew-Hamiltonian, and 

(4) G{s)^ = —G{s) skew-symmetric {n must be even). 

Symmetries (l)-(2) were studied in [12] and occur naturally in systems theory [HI E]- 
Stabilization of symmetric systems (1) by real symmetric output feedback was considered 
in [20], who showed that there may be no real feedback laws placing n real poles when 
n > m. The last two symmetries are natural to consider from the point of view of algebraic 
geometry. Theorem [3] gives an example of an m-input m-output skew-symmetric linear 
system of McMillan degree 2(™) such that every feedback law is real when placing real 
poles, showing that it is feasible to place real poles with real skew-symmetric feedback. 

Let A be a nondegenerate bilinear form on C^"^. The annihilator H^-'^ of an m-plane 
H in C^™ is also an m-plane. The association H i— )■ H^^ defines an involution la on 
the Grassmannian G{m,2m). If A is skew-symmetric, then the set of fixed points of la 
is the Lagrangian Grassmannian, LG{m), which is a manifold of dimension ('"^^). If A 
is symmetric, then the set of fixed points of la has two isomorphic components, either 
of which forms the (^) -dimensional orthogonal Grassmannian, OG{m), also called the 
spinor variety [TT] . 

It is classical (e.g., proved by Gaussian elimination [1]) that any two invertible complex 
symmetric matrices A, B are (transposae) congruent: there exists an invertible matrix X 
such that X"^ AX = B. Similarly, any two invertible complex skew-symmetric matrices 
are congruent. Thus, we will always assume that our forms are {x,y) = Ay, where A 
is either 



O 



2m 



Im 

Im. 



or J' 



2m 



Im 

-Im. 



in which Im is the m x m identity matrix. We omit the subscripts m and 2m when 
the dimensions are clear from context. A general linear subspace H G G{m, 2m) is the 
row space of a matrix of the form [Im '■ F], where F is an m x m matrix. A calculation 
shows that lj{II) is spanned by [J^ : F^^] and LoiH) is spanned by [Im '■ — F'^], so that 
H G LG{m) if and only if F is symmetric and H G OG{m) if and only if F is skew- 
symmetric. 

We reach the same conclusion for H of the form [F : Such isotropic planes form 
dense open subsets of LG{m) and OG{m). 

If we associate an m-input m-output transfer function G{s) of McMillan degree n to 
the row space of the matrix 

[Im : G{s)], 

we obtain a map 7: — > G(m, 2m) of degree n. The image is the Hermann-Martin 
curve [2T] of G{s). The set of all such transfer functions forms a dense open subset in the 
space of rational curves of degree n in the Grassmannian G{m,2m) [25]. Our first main 
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result identifies the sets of transfer functions with symmetries as natural subvarieties of 
the space of rational curves in the Grassmannian. 

Theorem 1. (1) The set of symmetric linear systems with m inputs and m outputs of 
McMillan degree n is a quasiprojective manifold of dimension (m+l)n + ('"^^) • It is a 
dense open subset of the space of rational curves of degree n in LG{m). 

(2) The set of Hamiltonian linear systems with m inputs and m outputs of even McMil- 
lan degree n is a quasiprojective manifold of dimension mn + ('"^^) • It is a subset of the 
space of rational curves 7 in the Grassmannian G{m, 2m) that satisfy: 

(1.1) 7(-^) = Ulis)). 

(3) The set of skew- Hamiltonian linear systems with m inputs and m outputs of McMil- 
lan degree n is a quasiprojective manifold of dimension mn + (™) . It is a subset of the 
space of rational curves 7 in the Grassmannian G{m, 2m) that satisfy: 

(1.2) 7(-s) = ^o{l{s)). 

(4) The set of skew- symmetric linear systems with m inputs and m outputs of even 
McMillan degree n = 2i is a quasiprojective manifold of dimension {m—l)n + (™) . It is 
naturally a dense open subset of the space of rational curves of degree C. in OG{m). 

The proof of Theorem [1] is straightforward and given in Section [21 following a proof of 
a version of the Kalman Realization Theorem p!B| Theorem 6.2-4] for symmetric transfer 
functions. We do not know if Hamiltonian systems form dense open subsets of the space 
of curves satisfying (11.11) . or if skew-Hamiltonian systems form dense open subsets of the 
space of curves satisfying (II. 2p . for these spaces of curves have yet to be studied. 

As sets of linear systems with symmetry are identified with quasiprojective algebraic 
varieties, the notion of genericity for complex systems makes sense. That is, a property 
is generic if it holds on a nonempty Zariski open subset of the corresponding space. 

Because symmetric and skew-symmetric transfer functions are open subsets of the mod- 
uli spaces of rational curves in the Lagrangian Grassmannian and orthogonal Grassman- 
nian, respectively, output feedback control by either static or dynamic symmetric and 
skew-symmetric linear systems is related to Schubert calculus on these Grassmannians, 
both classical (for static feedback laws) and quantum (for dynamic feedback). The main 
result of [12] concerned static symmetric feedback. We establish the analogous result for 
static skew-symmetric feedback. 

The symmetry of skew-Hamiltonian and skew-symmetric linear systems is preserved by 
static skew-symmetric output feedback, so it is natural to place poles with static skew- 
symmetric controllers. The poles of a skew-Hamiltonian transfer function are invariant 
under multiplication by —1, so there are essentially only [n/2\ poles to place. Similarly, 
poles of a skew- symmetric linear system occur with even multiplicity, and therefore a 
skew-symmetric linear system of even McMillan degree n has only n/2 poles to place. 
Our second main theorem gives necessary and sufficient conditions for pole placement 
with static skew-symmetric feedback. 
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Theorem 2. A general strictly proper skew- symmetric (respectively skew-Hamiltonian) 
transfer function G{s) withm inputs, m outputs, and McMillan degree n is pole- assignable 
with complex static skew- symmetric feedback compensators if and only if i = \ji/2\ < (™) . 

In particular, when m = 4, we see that a skew-symmetric system of McMillan degree 12 
or less is pole-assignable with complex static skew- symmetric feedback compensators, and 
skew-Hamiltonian systems of McMillan degree 13 or less are pole-assignable. Our third 
main result counts the number of feedback laws for a general skew-symmetric system of 
McMillan degree 2 . It also shows that there exist systems with real feedback laws, in 
a strong way — for these systems, every feedback law placing real poles is real. 

Theorem 3. For a general skew- symmetric linear system of McMillan degree 2(™), there 
will be exactly 

/mV 0!1!-- -(771-2)! 

" V2 ;'l!3!---(2m-3)! 

static complex skew- symmetric controllers that place a given general set of (™) poles. 

Moreover, for every m, there exists a skew-symmetric linear system with m inputs and 
outputs and McMillan degree 2 (^) such that for every choice of (™) real poles, there are 
are dm feedback laws, and every one is real. 

We prove these theorems in Section [3l The argument for Theorem [2] is influenced by 
the proof in [12] , but it is a considerable simplification. The skew-symmetric system with 
dm real feedback laws comes from the Wronski map in the Schubert calculus, and the 
result on reality is a restatement of a theorem of Purbhoo [23] . 

We do not address questions about dynamic feedback. If we use a dynamic compensator 
of McMillan degree q to place the poles of a linear system of McMillan degree n with one of 
these symmetries, then a calculation shows that the resulting system (of McMillan degree 
n -\- q) has the same symmetry as the original system if and only if the compensator 
had that same symmetry. Thus it is natural to consider dynamic control when both the 
system and compensator have the same symmetry. A dimension count gives the necessary 
condition that n + g be at most 

/m + l\ /m + l\ / m\ fm 

(m+l)g + I , mg + I , mg + I , and [m—\)q^ 



for generic pole placement of symmetric, Hamiltonian, skew-Hamiltonian, and skew- 
symmetric linear systems by dynamic controllers of the same symmetry. We do not know 
if these conditions are sufficient — this requires the generic surjectivity of the corresponding 
pole-placement map. 

If this dimension condition is necessary and sufficient, then the quantum Schubert 
calculus for LG{rn) and OG{m) [TTJ [18] may be used to count the number of dynamic 
compensators for symmetric or skew-symmetric systems (also [28] for symmetric compen- 
sators). Counting dynamic compensators for Hamiltonian and skew-Hamiltonian systems 
requires a deeper study of the corresponding spaces of curves. 

Our results and analysis also apply to discrete-time linear systems with these symme- 
tries, in the same way that continuous-time transfer functions are related to discrete-time 
transfer functions when there are no symmetries. 
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2. Geometry of state-space realization with symmetries 

We study the state-space realizations of transfer functions with symmetry and identify 
the spaces of such transfer functions as certain subvarieties of the moduli spaces of rational 
curves in Grassmannians. Portions of this material are classical or can be found in [12], 
but we include some proofs for completeness. We work entirely over the complex numbers. 
Write for the transpose of a matrix X and X^^ for [X^^)~^ = {X~^)~^. Recall that a 
square matrix X is symmetric if X~^ = X and skew-symmetric if X~^ = —X. 
Let J be the 2i x 2i matrix, 

" /" 



J 



-I 



where I is the i x i identity matrix. Note that = — J = J^^. A 2i x 2i matrix X is 
Hamiltonian if XJ is symmetric and skcw-Hamiltonian if XJ is skew-symmetric. 

Let m and n be positive integers, which we assume are fixed throughout. We write i for 
Lr?,/2j . Suppose that we have a time-dependent complex linear system with inputs u G C™, 
outputs y G C™, and McMillan degree n. This has a minimal state-space realization: 

^ X X = Ax + Bu 

^ ' y = Cx + Du, 

where A G C"^'^, B G C"^"^, C G C™^'', and D G C"*^'". 

We eliminate the internal states x by taking the Laplace transform of (12.11) and working 
in the frequency domain. Write x{s), y{s), and u{s) for the Laplace transforms of the 

time-dependent functions x{t), y{t), and u{t). Since -^x = sx, we have: 

sx = Ax + Bu and y = Cx + Du . 

We solve the first equation for x to get x = {si — A)~^Bu, and then substitute into the 
second to obtain 

y = [C{sl - Ay^B + D]u, 

which expresses the linear system in the frequency domain. Defining 

G{s) := C{sl -A)~^B + D, 

this expression becomes y = G{s)u, and we call G{s) the transfer function of this system, 
which is a matrix-valued rational function. The poles of the transfer function G{s) are 
the eigenvalues of A. The poles govern the long-term behavior of the system when u = 0. 
The transfer function G{s) is strictly proper if D = 0. 

Given a strictly proper linear system, a static linear feedback law is given by an m x m 
matrix F, where we set u = Fy + v. The resulting linear system is 

(2.2) X = {A + BFG)x + Bv y = Gx , 
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and its transfer function has poles at the roots of 

^{s) := det{sl -{A + BFC)). 

These poles govern the long-term behavior of the open loop system obtained from (12.21) 
by setting f = 0. A fundamental problem is: When is it possible to choose F to obtain a 
given choice of monic polynomial ip{s)7 A system is polc-assignablc if the map F i— )■ ip[s) 
is dominant (its image is a dense subset of the set of polynomials ^p{s)). Our main result 
concerns the pole-assignability of generic linear systems with skew-Hamiltonian and skew- 
symmetric symmetry. We first study the spaces of linear systems with symmetry. 

The complex general linear group GL{n) of invertible n x n matrices acts on the space 
of realizations (12. ip via 

X.{A,B,C,D) I — > {X-^AX, X-^B, CX,D) , 

where X G GL{n). Since 

CX{sI -X~^AXY^X~^B + D = C{sI-Ay^B + D, 

every realization in a GL(n)-orbit has the same transfer function. If we restrict this action 
to the dense open set of controllable and observable (hence minimal) state space realiza- 
tions, then the Kalman Realization Theorem [T6l Theorem 6.2-4] identifies the orbits 
with transfer functions and shows that GL{n) acts without fixed points. In particular, if 
{A, B, C, D) and (a, /3, 7, 5) are both minimal state space realizations of the same transfer 
function, then there is a unique X G GL{n) such that 

X.{A,B,C,D) = («,/3,7,5). 

Our first goal is to extend these classical facts to transfer functions with symmetries. 

Definition 4. A transfer function G{s) is symmetric, Hamiltonian, skew-Hamiltonian, or 
skew-symmetric if for all s G C we have, 

G{sy = G{s), G{sy = G{-s), G{sy = -G(-s), or G(s)^ = -G{s) , 

respectively. We will call a transfer function symmetric if it exhibits one of these four 
symmetries. Context will remove the resulting ambiguity. 

State-space realizations may also have symmetries. 

Definition 5. A realization (12. ip is symmetric if A is symmetric, B = G~^ , and D is 
symmetric. Symmetric realizations have symmetric transfer functions: 

G^{s) = B^{sl - A'^y^G^ + = G{sl - Ay^B + D = G{s) . 

A realization (12. ip with n even is Hamiltonian if A is Hamiltonian, JB = G~^ , and 
D is symmetric. Note that A^ = JAJ and B^ = GJ. Hamiltonian realizations have 
Hamiltonian transfer functions: 

G{-sy = B^{-sI -A'^y^G^ + D'^ = GJiysJLJ - JAjy^JB + D'^ 

= G{sI-Ay^B + D^ = G{s). 
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A realization (12. ip is skew-Hamiltonian if A is skew-symmetric, B = C^, and D is skew- 
symmetric. Skew-Hamiltonian realizations have skew-Hamiltonian transfer functions: 

= -C{-sI + A)-^B + D = G{s). 

Finally, a realization (12. ip with n even is skew-symmetric if A is skew-Hamiltonian, 
JB = , and D is skew-symmetric. In this case, A~^ = —JAJ and B^ = CJ. Skew- 
symmetric realizations have skew-symmetric transfer functions: 

-G{sy = -B^{sl - A^y^C^ - D'^ 

= -GJ{-sJIJ + JAJ)-\-JB) + D = G{s). 

Remark 6. In symmetric and Hamiltonian realizations, the matrix A has the same 
type (symmetric or Hamiltonian, respectively), while for skew-Hamiltonian and skew- 
symmetric realizations, the matrix A has the opposite type; it is skew-symmetric or skew- 
Hamiltonian, respectively. 

While there is no a priori reason that a symmetric transfer function would have an 
appropriately symmetric minimal state-space realization, that is indeed the case. The 
orthogonal group 0{n) is the subgroup of GL{n) consisting of matrices X with X"^ X = I, 
and when n is even, the symplcctic group Sp{n) is the subgroup of GL{n) consisting of 
matrices X with X"^ JX = J. We establish the analog of the Kalman Realization Theorem 
for transfer functions with symmetry. 

Proposition 7. A transfer function has one of the symmetry types — symmetric, Hamil- 
tonian, skew-Hamiltonian, or skew- symmetric — if and only if it has a minimal state-space 
realization having the corresponding symmetry type. 

Furthermore, if {A, B,G, D) and (a,/3,7,5) are two such minimal state-space real- 
izations of the same transfer function, then the unique matrix X G GL{n) such that 
X.{A, B,C, D) = (a,/3,7,5) is orthogonal for symmetric and skew-Hamiltonian transfer 
functions, and symplectic for Hamiltonian and skew- symmetric transfer functions. 

Following the proof for symmetric transfer functions in [12] (see also [8]), we give the 
proof in the cases of skew-Hamiltonian and skew-symmetric transfer functions. The case 
of Hamiltonian transfer functions is similar, and is due to Brockett and Rahimi [1]. Also, 
the first half, concerning symmetric realizations, is due to Brockett [2]. 

Proof. We prove the forward implication in the first statement as we have already shown 
that a state-space realization having one of these symmetries gives a transfer function 
with the same symmetry. 

Suppose that G{s) = —G~^{—s) is a skew-Hamiltonian transfer function with minimal 
state-space realization {A, B,G, D). Since 

-G{-sy = -B^{-sI - A^y^C^ - D'^ = B^{sl + A'^y^G^ - D'^ , 
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{—A^, , , —D^) is also a minimal realization. By the Kalman Realization Theorem, 
there is a unique invertible matrix X such that 

(A 5, C, D) = X.{-A^, C^, B^, -D^) 

= {-X-^A^X, X-^C^, B^X, -D^) 

= {-X-\-X-'A'^XyX, X-\B'^XY, {X-^C'^yX, -D^) 

= {X-^X'^ AX-^X, X-^X^B, CX-^X, -D'^). 

= {X-'X^).iA, B, C, D). 

Since {A, B, C, D) = I. {A, B, C, D), another application of the Kalman Realization 
Theorem gives us that X~^X^ = /; thus, X is symmetric. 

An invertible complex symmetric matrix X admits a Tagaki factorization X = Y^Y, 
with Y invertible [14., Corollary 4.4.4]. Then the realization {YAY-^,YB,CY-^) of the 
transfer function G{s) is skew-Hamiltonian. Indeed, 

{YAY-^ = Y-^A'^Y^ = -Y-^X^AX-^Y^ 

= -Y-^Y^YAY-'^Y-^Y^ = -YAY~\ 
Similarly, (YB^ = CY-\ Indeed, C = B^X so B^ = CX-\ and thus 

{YBy = B^Y^ = CX-^Y^ = C{Y-^Y-^)Y^ = CY'^ . 

Now suppose that (A, 5, C, Z^) and (a,/3,7,(5) are minimal skew-Hamiltonian realiza- 
tions of the same transfer function. Let X G GLin) be the unique matrix such that 

X.{A,B,C,D) = {X'^AX,X-^B,CX,D) = (a,/3,7,5). 

Then we have 

a = -a^ = -X^A'^X-^ = X^ AX'^ , 

and similarly f3 = X^B and 7 = CX'^ so that X'^ .{A, B,C, D) = {a,f3,-f,6). It 
follows that X^^ = X^ , by the uniqueness of X. But then X~^X = I and so X G 0{n) 
is orthogonal. 

Consider next the case that G{s) = —G~^{s) is a skew-symmetric transfer function with 
minimal state-space realization {A, B,C, D). Since = , we see that —G{sy equals 

-B^{sl - A'^y^C^ - D'^ = B^J^J{-sI + A^)-'jJ^C^-D'^ 

= {JBY{-sJ-HJ-^ + JTaTjT)-1(C7J)T _ 

= {jbY{si -{-jajY)-\cjY - D'^ , 

and so {—{JAJY-, {CJY-, {JBY ■, ~D^) is also a minimal realization of G{s). Then there 
is a unique X G GL{n) such that 

(A, 5, C, D) = X.i-{JAJY, (cjY, (JbY, -d'^) 

= {-x-^{jajYx, x-\cjY, (JbYx, -d^). 

Substituting the equality into itself and simplifying, we see that 

(2.3) {A,B,C,D) = {X-^J^X^JAJX-^J^X, X'^J^X^JB, CJX'^J^X, D). 
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The right-hand side of (12. 3p is not immediately seen to have the form R.{A, B, C, D), but 
if we set R := JX~~^ X and use that — J"*" = J~~^ and = —J, we obtain 

= x-^j-^x^j-^ = x-^j^x^j, 

which shows that (Q is {R^^AR,R-^B,CR,D). We conclude that i? = J, so that 
-I = JX-^rX = J-^X-^JX, and so {JX^ = -JX is skew-symmetric. 

Paralleling the argument for skew-Hamiltonian transfer functions, we will use this to 
obtain a skew-symmetric realization. The key ingredient is a factorization of the matrix 
JX. An invertible skew-symmetric matrix Z admits a Tagaki-like factorization Z = 
JY. Indeed, Z = USU~^ for a unitary matrix U and a block diagonal S which is a 
direct sum of 2 x 2 blocks of the form [ -?a o ] with a e C \ {0} (e.g., see Problem 26 in 

Chapter 4.4]). Thus, after block scaling with blocks ^ ^ and applying a permutation 

similarity, we arrive at the claimed factorization. 

In the factorization JX = Y^ JY of the invertible skew-symmetric matrix JX, the 
matrix Y is also invertible. Then the realization {YAY~^,YB,CY~^,D) of the transfer 
function G{s) is skew-symmetric. Indeed, as D = —D^ and we have J = —J~^ = J~~^ , 
X^J^ = -JX = -Y^JY, and J-^X'^ = Y'^JY'^ , we obtain 

{YAY-^jy = J^Y-^A'^Y^ = -~j'^Y-^{X-\jAjyX)'^Y^ 

= -J^Y-^X^ JAJX-^Y^ = J'^Y-^ X^ J'^ AJ-^ X-^Y^ 
= -J^Y-^Y^ JYAY-^JY-^Y^ = -YAY'^J. 

Similarly, CY'^ = {JB^XY'^ = (JYB)^, so the realization {YAY-\YB,CY-\D) of 
G{s) is skew-symmetric. 

Finally, suppose that {A, B,C, D) and (a,/3,7,5) are minimal skew-symmetric realiza- 
tions of the same transfer function. Let X G GL{n) be the unique matrix such that 

X.{A,B,C,D) = {X-^AX,X-^B,CX,D) = (a,/3,7,5). 

Recall that — = JAJ and a = ~Ja^ J, and so 

a = -Ja^J = -JX^A'^X-^J = JX^ JAJX'^ J = J-^X^ JAJ-^X'^ J . 

Similarly recall that = JB and /3 = J7^, so that 

P = JT^T _ j^[CXy = J'^X^C^ = -JX^JB = J-^X^JB. 

Since we also have 7 = CJ~^X^^ J, we see that (a,/?, 7, 5) = R.{A, B,C, D), where 
R = J~^X~^ J. By the uniqueness of X, we have X = J~^X^^ J so that X^ JX = J and 
so X G Sp{n) is symplectic. □ 

We use this proposition to compute the dimensions of the corresponding spaces of 
transfer functions/rational curves, the first part of the proof of Theorem [TJ 

Corollary 8. The set of transfer functions with a fixed symmetry is a smooth quasipro- 
jective algebraic manifold. For symmetric, Hamiltonian, skew-Hamiltonian, and skew- 
symmetric transfer functions, these have respective dimensions 

^m+l)n + (™+') , mn + ("+^) , mn + (^) , and {m-l)n + (™) . 
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Proof. The space of minimal symmetric realizations is a Zariski-open subset of an affine 
space. By Proposition [TJ the set of transfer functions with a fixed symmetry type is 
identified with the set of orbits of an algebraic group (0(n) or Sp{n)) acting freely on the 
space of minimal symmetric realizations, which is an open subset of a vector space. The 
first statement of the corollary follows from this as the set of such orbits has a natural 
structure as a smooth algebraic variety [ISl Th. 9.16]. 

For the second, we note that the dimension of the orbit space is the difference of the 
dimensions of the space of symmetric realizations and of the group. The orthogonal group 
0{n) has dimension (2) and the symplectic group Sp{n) has dimension ("^^) [ID]. Spaces 
of symmetric and Hamiltonian realizations both have the same dimension 

fn+l\ fm+l\ 

while the spaces of skew-Hamiltonian and skew-symmetric realizations both have dimen- 
sion 

f n\ ( m\ 

-|- nm -|- 

The corollary now follows. □ 

We now complete the proof of Theorem [T] by identifying the Hermann-Martin curves of 
symmetric and skew-symmetric transfer functions with dense open subsets of the appro- 
priate spaces of rational curves in the Lagrangian and Orthogonal Grassmannians. First 
note that if G{s) is symmetric (respectively skew-symmetric) then for s G P^, the row 
space K{s) of the matrix [/„ : G{s)] lies in LG{m) (respectively in OG{m)). Thus the 
Hermann-Martin curve is a curve in LG{m) (respectively in OG{m)). 

To finish, we show that that these Hermann-Martin curves are dense in the correspond- 
ing spaces of rational curves, which is a consequence of their having the same dimension. 
The space of rational curves in LG{m) of degree d has dimension d{m+l) + ('"^^) and 
the space of rational curves in OG{m) of degree d has dimension 2d{m—l) + [T71 [T5] . 
Thus Theorem [T] follows if we knew that a curve in LG{m) of McMillan degree n has 
degree n in LG{m) and a curve in OG{m) of McMillan degree n = 2i has degree i in 
OG{m). 

These facts are well-known. The inclusion LG{m) )■ G{m, 2m) arises from a linear 
map on the standard projective embedding of LG{m), so rational curves of degree n in 
LG{m) have degree n in G{m,2m). However, the inclusion OG{m) ^ G{m,2m) arises 
from the second Veronese map on the natural projective embedding of OG{m). Thus a 
rational curve of degree i in OG{m) will have degree 2i in G{m, 2m), and this completes 
the proof of Theorem [H 

3. Static skew-symmetric state feedback control 

Suppose that we have a strictly proper {D = 0) skew-Hamiltonian or skew-symmetric 
linear system with a minimal state-space realization 



(3.1) 



X = Ax + Bu y = Cx . 
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If we introduce a static linear state-space feedback law u = Fy + v, then the new system 

(3.2) X = {A + BFC)x + Bv y = Cx 

has the same symmetry as the original system when F is skew-symmetric. (The elemen- 
tary calculation is given below.) We investigate the control of such linear systems with 
complex skew-symmetric static state-space feedback. We first establish the necessary and 
sufficient conditions for generic pole placement of Theorem [21 relate skew-symmetric feed- 
back control to the Schubert calculus on the orthogonal Grassmannian, and then prove 
Theorem 121 counting the number of feedback laws that place a generic set of poles of 
a generic skew-symmetric transfer function with McMillan degree 2(^). We do not yet 
know how to count the controllers of a skew-Hamiltonian transfer function of McMillan 
degree n when (™) = [n/2\. 

Proof of Theorem [H We give the proof for generic pole-assignability of skew-symmetric 
transfer functions and indicate how the argument changes for skew-Hamiltonian transfer 
functions. This follows and simplifies the arguments in [T^ . 

We identify skew-symmetric N x N matrices with the vector space A^C'^, where the 
elementary decomposable tensor Cj A Cj {i 7^ j) corresponds to the matrix having 1 in 
position —1 in position and in other positions. 

Suppose that (13.11) is skew-symmetric, so that A is skew-Hamiltonian, {AJ)~^ = —AJ, 
and C = J. If we have a feedback law u = Fy + v where F is skew-symmetric, then 
the new system (13. 2p remains skew-symmetric as A + BFB^ J is skew-Hamiltonian. Thus, 
the characteristic polynomial 

(3.3) ip{s) := det{sl - {A + BFB^J)) = det{sJ-{AJ-BFB^)), 

is the determinant of a skew-symmetric matrix and is therefore a square (its determinant 
is the square of its Pfaffian). Thus it is natural to ask for skew-symmetric feedback laws 
F which place these i = n/2 roots (which are poles of the transfer function). 

The pole placement map sends a skew-symmetric matrix F G A^C™' to the degree 2i 
polynomial v^(s). Since this monic polynomial is a square, its last i coefficients (those 
of s^^~^, . . . , s^) determine its first i coefficients. These coefficients are, up to a sign, 
the elementary symmetric functions of the eigenvalues of A + BFB^ J. By the Newton 
identities, these coefficients determine, and are determined by, the Newton power sums 
which are the traces of {A + BFB^ JY for A; = 1, . . . , 

To show generic pole-assignability, we only need to exhibit one choice of matrices A, B 
for which the map 

^ : A^ 3 F ^ (Tr((A + BFB'^ jf) | = 1, . . . , £) G 

is dominant. We do this by showing that the differential d^Q at G A^C™ is surjective. 

Let be distinct numbers and . . . , /3m be numbers such that the (™) prod- 

ucts I3i(3j for i < j are distinct, and such that /3f 7^ for every i 7^ j. Let D = 
diag(ai, . . . , ae) be the diagonal matrix with entries ai, . . . ,ae, and let A be the block 
diagonal matrix [{j* j^]. Finally, let B be the matrix with entries for i = 1, . . . ,2i and 
j = 1, . . . , m. For this choice of A, B, the differential d"^ is surjective at G A^C". 
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To see surjectivity, note that the differential at is the linear map 

(3.4) d'^o ■■ F ^ {k- Tt{A''-^ ■ BFB^J) \k = !,...,£). 

Consider this map on the basis element Cj A Cj of A^C™. A direct calculation shows that 
B{ei A ej)B^ is the vector hi A 6j, where bi, . . . ,bm are the columns of B. For our choice 
of B, the {p, g)-entry of hi A bj is 

(A/3,)^-'(/3r^-/3r') 



det 

Then the map d^Q sends the vector Cj A Cj to the vector 

21 

'kY,{A''~\h^h)J)p,p\k = l,... 
p=i 

Since {{h A bj)J)^^^ = {{h A bj)j) ^^^^^^^ = {h A bj)p+e,p, this vector is 

'2i53aJ-i(ftft)''-'(/?,'-,3j 



p=l 



Thus, (i\l/o is represented by the i x (™) matrix which is the product of the matrices 

diag(2, 4, . . . , 2£) ■ {a'p-Xzl::', ■ {mrTr<l<j'L ■ diag(/5f - < ^ < J < m) , 

and so its rank is the minimum of £ and (™), which proves Theorem [2] for skew-symmetric 
linear systems. 

Suppose now that (13. ip is skew-Hamiltonian, so that A is skew-symmetric and C = B^ . 
Under a skew-symmetric feedback law u = Fy + v, the new system (13. 2p remains skew- 
Hamiltonian as A + BFB^ is skew-symmetric. The characteristic polynomial 

(3.5) ¥?(s) := det(sJ - (A + SFB^)) 

satisfies ^{s) = (— l)"(y9(— s), so its nonzero roots A occur in pairs ±A. Thus it is natural 
to ask for skew-symmetric feedback laws F which place these i = \n/2\ pairs of roots 
(which are poles of the transfer function). 

The pole placement map sends a skew-symmetric matrix F G A^C™ to the degree n 
polynomial ip{s). As before, we investigate the surjectivity of the pole-placement map by 
considering the map associating Newton power sums, which is 

^ : A^ 9 F I — y {Tt{A - BFB'^f \ k = l,...,n) G 

Since A + BFB^ J is skew-symmetric, the trace is zero if k is odd, which is why the 
codomain of this map is C^. Thus, the differential of at is 

d^o : F I — > ~{2k ■ Ti{A'^''-^BFB^) \k = l,...J). 

Let D be the same diagonal matrix as before. If n is even, let A be the block matrix [ _£) o* ] 
and B be the same as before. If n is odd, then add a new first row and column of Os to A 
and extend B with the row {Pf^ \ i = 1, . . . ,m). Then nearly the same calculation as before 
shows that d'^o is surjective when i < which completes the proof of Theorem [2l □ 
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Before proving Theorem [3l we first recall some standard matrix manipulations which 
transform the problem of finding matrices F which place the poles of the transfer func- 
tion (13 ■2p into a geometric problem on a Grassmannian. We then make some definitions. 

These poles are the roots of the characteristic polynomial ^{s) = det(sJ„ — {A + BFC)) 
of the matrix in (13. 2p . The rational function f{s)/ det(s/„ — A) equals the determinant 
of the product 





{sin -A)-' 





" 




sin 


-A- 


-BFC 


BF 


-B' 


'In 


" 






-C{sln-A)-' 


Im 
















Im 





C Im 















/ 


m 


















Im 


F 


Im 




which is 


































{sin -A)-' 


o" 




sin 


-A 





-B 






'In 





-{sin - 




'B 




-C{sln-A)-^ 
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Im 












Im 


C{sln - 




'B 









Im 






F 





Im 









F 


Im 





The transfer function G{s) = C{sln — A) ^B admits a left coprime factorization into 
matrices of polynomials D{s)~^N{s), where det D{s) = det(s/„ — A), and so we have 



(3.6) 



ip{s) = det(s/„ — A) det 



Im 

F 



G{s) 



det 



D{s) 
F 



N{s) 



Thus, s is a pole of the transfer function of (13. 2p if and only if the matrix on the right 
of (13.61) does not have full rank. Geometrically, if K{s) is the row space of [D{s) : N{s)] 
and H is the row space of [F : which are both m-planes in C^"*, then K{s)r\II ^ {0}. 

It follows that the feedback laws F which place a given set of poles si, . . . , s„ correspond 
to those H G OG{m) such that 



(3.7) 



K{si) nH 7^ {0} , for each i = 1, 



n. 



When G{s) is skew-symmetric, then n = 2i and ip{s) has only i distinct roots, say 
Si, . . . ,Si. In this case, we also have that K{si) is isotropic and the set of if G OG{m) 
satisfying (13.71) defines a Schubert subvariety of OG{m). A consequence of the surjectivity 
of the map d^Q (13.41) for general A, B is that these Schubert varieties meet generically 
transversally (in the open set consisting of H of the form [F : /m])- Thus, when ^ = (™) 
there are finitely many H satisfying (13. 7p for i = !,...,£, and their count is bounded 
above by an intersection number computed in the Schubert calculus on OG{m) [H] to be 



2 j'l!3! 



0!1! ■ ■ - (m - 2)! 



(2m -3)! 

We complete the proof of Theorem [3] by exhibiting a specific skew-symmetric transfer 
function G{s) of McMillan degree 2(™) such that there are exactly skew-symmetric 
feedback laws placing any given ('^) real poles, and all the feedback laws are real. This 
also shows that generic systems have dm feedback laws. 

The argument uses a result of Purbhoo [23] concerning the reality of the Wronski map, 
which we will transfer into the language of systems theory. Suppose that C^™ has ordered 
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basis ei, . . . , e2ri 



Let 7(5) be the vector- valued function 7: C — t- C^™ with components 

2 ^m-l ^ f-s)2m-2 /^Nm (_s)"^^^ 1 



(m-1)! 72 ' (2m - 2)! ' ■ ■ ■ 

for z = 1, . . . , m— 1 and fm(s 



)• 



If we set v,{s) := (^)' 7(5, 

(— l)™~^e2m), then the row span K{s) of 

m = 5, giving a 5 x 10 matrix with rows Vi{s 



(m-1)! 

(f)'""SW + ;^(«.» + 

fm(s) is isotropic. We display this for 
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V2 4! 
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1 s3 
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7! 


6! 


5! 


4! 
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3! 
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V2 


4! 


s3 
3! 


2! 


— S 






If we write K(s) = [D(s) : A^(s)], then D{s) is upper triangular with diagonal (1, . . . , 1, -\/2), 



_2£_ 
5 7! 

2 

5 6! 

3 5! 
V2 4! 



Set G'(s) 
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»-^iv(.; 
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1 
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Theorem |3] follows from the following facts about the transfer function G{s). 

Proposition 9. The transfer function G{s) is skew- symmetric with McMillan degree 
2(™). Any set of (™) distinct real poles is placed by exactly dm skew- symmetric feedback 
laws, with each one real. 

Proof. Since the isotropic m-plane K{s) is the the row space of [/„ : G{s)], we conclude 
that G{s) is skew-symmetric. Let V ~ (CSm-i ^ (^2m ^^le subspace with ordered basis 



ei, 



-1, e„ 



-l)"-^e2. 



2,e 



7 "= m+l; 



e2m-l. 



The nondegenerate symmetric bilinear form on C^"* restricts to a nondegenerate symmet- 
ric bilinear form on V, and the map H W := H (IV sends a maximal isotropic subspace 
H of C^™ to a maximal isotropic subspace of V, inducing an isomorphism between OG{m) 
and the space BOG{m—l) of maximal isotropic subspaces of V c:^ C^™""^. The reason for 
this is that for each W G BOG{m—l) there are two maximal isotropic subspaces H of 
£2m. containing W, exactly one of which lies in OG{m). When W is real, both isotropic 
subspaces H containing W are also real. 

Also, 7(5) is a rational normal curve in V , as it involves the monomials 1, . . . , s^*""^. 
Furthermore, L(s) := K{s) fl is the (m— l)-plane osculating 7(5), and L(s) is isotropic. 
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The problem of studying isotropic subspaces W of V that meet r = (^) osculating 
planes L{si), . . . , L{sr) was studied by Purbhoo [23] in the context of the Wronski map 
from BOG{m—l) = OG{m), which extends the pole placement map from [F : Im] ^ (p{s) 
for the transfer function G{s). This map is surjective onto the space of polynomials of 
degree 2r which are squares of polynomials, and it has finite fibers of algebraic degree dm- 
This implies that there are at most dm feedback laws placing a given set of r poles. It 
also implies that any given isotropic plane H meets at most r subspaces of the form L{s), 
including L(oo) = [0 : Im-i]- 

Purbhoo [221 Theorem 3] showed that if Si, . . . , were real, then there are exactly dm 
real isotropic planes W in BOG{m — 1) such that L{si) nW {0}, for each i = 1, . . . ,r. 
For each such W, let H be the unique isotropic plane in OG{m) containing W, which is 
necessarily real. Then H r\ K{si) ^ for each z, and so H corresponds to a real feedback 
law if H has the form \lm '■ F\. But this is guaranteed for otherwise HnK{oo) ^ 0, which 
would imply W fl L(oo) {0}, which is impossible as H already meets the maximum 
number of subspaces of the form L{s). Lastly, the transfer function has McMillan degree 
2r since the image of the pole placement map (a linear projection) meets the set of 
polynomials of this degree. □ 
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